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ABSTRACT 

The nucleus of the bulge-dominated, multiply-barred SO/a galaxy NGC 2681 is stud- 
ied in detail, using high resolution Hubble Space Telescope FOC and NICMOS imaging 
and FOS spectroscopy. The ionised gas central velocity dispersion is found to increase 
by a factor ~ 2 when narrowing the aperture from R ~ 1'.'5 (ground) to R ~ Of.'l (FOS). 
Dynamical modeling of these velocity dispersions suggests that NGC 2681 does host 
a supermassive black hole (BH) for which one can estimate a firm mass upper limit 
M. < 6 x 10 7 M o . This upper limit is consistent with the relation between the central 
BH mass and velocity dispersion M, — a known for other galaxies. The emission line 
ratios place the nucleus of NGC 2681 among LINERs. It is likely that the emission line 
region comes from a rather mild, but steady, feeding of gas to the central BH in this 
galaxy. The inner stellar population lacks any measurable color gradient (to a radius of 
0.6 kpc) from the infrared to the ultraviolet, consistently with FOC, FOS and IUE data, 
all indicating that this system underwent a starburst ~ 1 Gyr ago that encompassed 
its whole interior, down to its very center. The most likely source of such a widely- 
distributed starburst is the dumping of tidally-extruded gas from a galaxy neighbor. 
If so, then NGC 2681 can be considered as the older brother of M82, seen face-on as 
opposed to the edge-on view we have for M82. 

Subject headings: galaxies: individual (NGC 2681) — galaxies: nuclei — galaxies: 
photometry — galaxies: spiral 
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1. INTRODUCTION 

Hubble Space Telescope (HST), Faint Object Camera (FOC) images for the central regions 
of the SO/a galaxy NGC 2681 were obtained in four ultraviolet bands in 1993, in the context of 
a project aimed at understanding the origin of UV emission from bulge-dominated galaxies. A 
point-like source — which we term spike for short — was evident at its center with a photometric 
profile being indistinguishable from the PSF of the aberrated, pre-COSTAR HST images (Bertola 
et al. 1995). 

Motivated by the related discovery of the variable nature of the UV-bright spike in NGC 4552, 
whose flare has been interpreted as an accretion event onto a central supermassive black hole (BH) 
(Renzini et al. 1995; Cappellari et al. 1999), we obtained additional FOC imaging and Faint Object 
Spectrograph (FOS) spectra of the NGC 2681 spike in 1997, this time with the COSTAR-corrected 
optics of HST. Unlike the case of NGC 4552, this latter set of HST observations did not reveal any 
variability of the central spike in NGC 2681 (not surprising, as we show here that a possible UV 
flare of the same luminosity as seen in NGC 4552 would have escaped detection in NGC 2681 due 
to its steeper surface brightness profile). 

In this paper we demonstrate that for this galaxy convincing evidence for a central BH comes 
instead from the modeling of the circumnuclear ionised gas kinematics derived from the FOS spec- 
trum, properly combined with FOC photometry. In so doing, we find we can also place interesting 
constraints on origin of the young stellar population in the center of this galaxy. Section 2 presents 
the known information about NGC 2681 and the observations we obtained. Section 3 presents the 
spectroscopic observations. Section 4 discusses the analysis of photometric HST data. Section 5 
derives the physical parameters for the galaxy inner nucleus, including our best estimate for upper 
limit to the mass of the central BH that resides in this galaxy. Section 6 discusses and summarizes 
our results. 

2. OBSERVATIONS AND DATA 

2.1. Existing Data for NGC 2681 

NGC 2681 is a bulge-dominated SO/a galaxy at an estimated distance of 13.3 Mpc (Tully 1988), 
corresponding to a linear scale of ~ 64 pc arcsec^ 1 and an absolute B magnitude Mb — —19.74 
(de Vaucouleurs et al. 1991). The value of the central magnesium index (Mg2=0.15; Burstein et al. 
1988; Jablonka et al. 1996) is low compared to values typical of spheroids of comparable luminosity. 
The IUE spectrum of this galaxy (Burstein et al. 1988) shows that the stellar population in its 
central region has a ultraviolet spectrum slowly decreasing towards shorter wavelengths, which is 
interpreted as evidence for a recent 1 Gyr) burst of star formation (Burstein et al. 1988). 

NGC 2681 has two concentric, misaligned bars (Wozniak et al. 1995) as seen in ground based 
images. A third inner bar is seen in this galaxy in HST/NICMOS images (Erwin & Sparke 1999) 
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. This innermost bar is also seen in our higher resolution FOC/96 images. Note that NGC 2681 is 
viewed almost face-on and thus these structures cannot be due to projection effects (i.e. we cannot 
misinterpret a change of flattening in an axisymmetric structure for a change of ellipticity due to 
a bar). 

Ground based B — I color maps do not show evidence of a ring or a grand design spiral, 
although two barely visible, small spiral arms exist in the disk. Based on our FOC/96 images the, 
spiral structure pointed out by Wozniak et al. (1995) appears to extend to within the central few 
arcsecs of this galaxy. Galactic absorption for this galaxy is E(B — V) ~ 0.02 (Burstein & Heiles 
1984; Schlegel, Finkbeiner, & Davis 1998). Ground-based Ha images also reveal a few H II regions 
and strong nuclear emission (Wozniak et al. 1995). Boker et al. (1999) has reported detection of 
Centrally-concentrated emission (r < 2"5) via Pa a imaging using HST/NICMOS. On the basis 
of their spectroscopic optical observations, Ho et al. (1997) classify the nucleus of this galaxy as a 
LINER. 

2.2. HST Observations 

FOS/red detector spectra of the nucleus of NGC 2681 were obtained in 1997, using the / .'21 x 
/ /21 square aperture (0.25-PAIR). A complete log of FOS observations is given in Table 1. Analysis 
of the FOS acquisition steps assures us that the nucleus of NGC 2681 was correctly centered in the 
FOS aperture (with a 0"025 accuracy). This is confirmed by the fact that the central flux measured 
within the FOS aperture on our COSTAR-Corrected F342W FOC/96 image agrees with the flux 
predicted for that FOC filter from the FOS G270H spectrum. 

FOC observations of NGC 2681 include our own 1993 Pre-COSTAR and subsequent 1997 
COSTAR-corrected, FOC/96 images. A detailed summary of the FOC instrumental configurations 
and exposure times is given in Table 2. In this work we also make use of an archival NICMOS image 
taken on 1997 June 7 with the F160W (space-based H-band), during one of the NIC3 "campaigns" 
(Boker et al. 1999). The scale of this image is C//2 pixel -1 and the field of view (FoV) is 51" x 51". 
All of the spectra and images used in this work have been de-archived and recalibrated using the 
most recent version of the STScI data pipeline and the latest calibration files. 

3. FOS SPECTROSCOPY 

3.1. The Overall Features of the FOS Spectrum 

The combined FOS nuclear spectrum of NGC 2681 is presented in Figure 1. Clearly recog- 
nizable in absorption are the Mg II A2800, Mg I A2852, and the Balmer lines H5, H^y and Hf3, as 
well as Ha in emission. Strong Balmer absorption lines in the spectrum of a galaxy are, of course, 
indicative of a relatively young stellar component (1-2 Gyr). This is consistent with the relative 
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absorption line strengths of Mg II and Mg I, which are indicative of a late A-type main sequence 
for the stars producing this spectrum (Wu et al. 1983), as well as with the interpretation of IUE 
spectrum of this galaxy (Burstein et al. 1988). Also prominent are emission lines of the C II] A2326 
multiplet, [0 II] A3727, [0 I] A6300, [N II] AA6548,6584, and [S II] AA6716,6731. 

Figure 2 shows a comparison between the nuclear FOS spectrum of NGC 2681 and the spectrum 
for this galaxy within the 10" x 20" IUE aperture, both normalized to their respective visual fluxes. 
The IUE aperture covers an area pa 4000 times larger than the FOS aperture we used, and the FOS 
flux contributes pa 4% of the IUE flux. This can be verified by comparing the UV flux in the FOS 
and IUE spectra (Burstein et al. 1988). Yet, the two spectra define the same kind of relatively 
young stellar population within combined errors of the IUE spectrum (0.25 mag maximum error 
quoted by Burstein et al.) and the 5% relative error of the FOS spectrum as quoted in the latest 
HST Data Handbook (Version 3). 

The luminosity-weighted average radius within the FOS aperture is R a ~ / /06, while for the 
IUE aperture it is R a pa 5"2. The agreement between the IUE and FOS spectra thus gives strong 
evidence that little, or no color gradient exists for the stellar population in the center region of 
NGC 2681. As a result, we conclude that whatever redleak does exist for the UV FOC images, the 
lack of a color gradient justifies fitting the same two-dimensional light distribution to all the FOC 
images and the NICMOS image (Section 4). Interestingly, both the lack of any perceptible dust 
lanes seen in the FOC images, combined with the lack of a color gradient seen in the IUE to FOS 
comparison, strongly suggest that the inner 0.6 kpc of NGC 2681 contains little dust. 

3.2. The Emission Lines 

The emission lines in our spectra have been analyzed by a "template subtraction" procedure 
as in Ho et al. (1997). Starting from the high S/N spectra of Ho et al. (1995), we have used 
a minimization algorithm to construct an emission-line- free template for their NGC 2681 2"x4" 
aperture spectrum. The best match was obtained using a combination of their spectra of NGC 205 
and NGC 4339. We use this template to subtract the underlying stellar contribution from the lower 
S/N 0"2x0 / /2 FOS spectrum. The upper panel of Figure 3 illustrates the template subtraction from 
the G650L spectrum (lacking any reliable emission), while the lower panel shows both the template 
subtraction and the subsequent fit to the prominent emission lines of the template-subtracted FOS 
G780H spectrum. 

The emission lines in the continuum-subtracted optical FOS spectrum have been modeled using 
the Levenberg-Marquardt algorithm (Press et al. 1992) to fit a sum of Gaussian components to the 
observed spectrum. We adopt the minimum set of free parameters able to generate a satisfactory 
fit by assuming that all emission lines can be fit by a Gaussian profile with the same redshift and 
the same velocity width within a given grating. The derived parameters for the emission lines are 
given in Table 3. The la formal errors of the line parameters are estimated from the covariance 
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matrix. The C II] emission line has also been fit in the G270H spectrum (Table 3) by assuming a 
linear template under the emission feature. 

The integrated Ha+[N II] flux in our C//21 x (X'21 aperture is only a factor ps 1.9 smaller than 
the flux measured by Ho et al. (1997) within a 2" x 4". Given that the ratio of the two apertures 
is a factor ~ 180, this small difference in the measured flux shows that the gaseous emission is 
very centrally concentrated and, in particular, comes mainly from a region not exceeding in size 
the FOS aperture. 



3.3. Emission Line Diagnostics 

Figure 4 compares the emission line ratios obtained from our HST spectra of the nuclei of 
NGC 2681 and NGC 4552 (Cappellari et al. 1999), to those for Seyfert galaxies, LINERs, and 
H II galaxies from Ho et al. (1997). Note that the [O III] and H/3 fluxes were not measured from 
our G650L FOS spectrum, as these lines have low S/N in our data. Rather, for these emission 
lines, we resort to the Ho et al. (1997) ground-based values. This choice cannot affect the physical 
classification of NGC 2681 nucleus as most of the flux observed by Ho et al. comes from a region 
within our HST aperture and — what is more — the basic discrimination on the AGN/non-AGN 
nature of the nucleus comes from our own (FOS) emission line ratios. 

As can be seen, the nucleus of NGC 2681 is clearly classified as a LINER from our FOS 
spectroscopy, in accordance with the ground based classification of Ho et al. (1997). We also note 
with interest that, despite the structural differences of the central regions of the two galaxies, the 
location of the innermost source of NGC 2681 is surprisingly close to that of the variable spike in 
NGC 4552 in these spectroscopic diagnostic diagrams. 



3.4. Constraints on the Size of the Emitting Region in NGC 2681 

A rough estimate of the size of the line-emitting region can be obtained from the density as 
set by the forbidden lines, and by the Ha luminosity (Osterbrock 1989): 

V ~ L(Ha) 

where V is the volume occupied by the emitting gas, L(Ha) is the total luminosity of the Ha line, 
/ is the volume filling factor, n e is the electron density, a^ a = 1.17 x 10~ 13 cm 3 s _1 , h is the Planck 
constant and vn a the frequency of the Ha line. 

Adopting n c ~ 1860 from Ho et al. (1997), a filling factor / = 1, a distance d = 13.3 Mpc 
(Tully 1988), and spherical geometry, we obtain L(Ha) « 2x 10 38 erg s" 1 and a size of the emitting 
region r « 1.2 pc (« / /02). This estimate is within the upper limit imposed by the size of our FOS 
aperture (« 6 pc). The latter size would correspond to a much lower filling factor of « 5 x 10~ 3 . 



-6- 



4. FOC AND NICMOS PHOTOMETRY 

Basic ingredients of the dynamical modeling we present later are the availability of a high 
resolution surface brightness profile of the galaxy, as well as observational support to the adopted 
hypothesis of a constant stellar mass-to- light (M/L) ratio. We show here that the data provided 
by the FOC and NICMOS images do assure both pieces of needed information. 

Ideally, in order to obtain reliable surface brightness profiles and color gradients, one should 
correct the FOC/96 images for nonlinearity effects in the FOC (Nota et al. 1996), deconvolve all im- 
ages with the PSF and, finally, derive the photometric profiles and the corresponding color gradients 
from the corrected frames. However, in our case this process is numerically very unstable, owing 
both to the very low signal-to-noise ratio of the FOC images and the approximations introduced 
by the nonlinearity correction which are greatly amplified by the deconvolution process. 

The similarity of the available FOS and IUE spectroscopy (Section 3) strongly suggests a lack 
of detectable color gradient between R rs 0'/06 and R « 5'/2. Moreover, the good agreement between 
the U (FOC/F342W) and H (NICMOS/F160W) profiles in the radial range 1"-10"— not affected 
by dark subtraction, nonlinearity and redleak effects — confirms the absence of any significant stellar 
population gradient in the central regions of this galaxy (see Figures 5, 6). For these reasons we 
feel permitted to construct a unique bi-dimensional photometric model of the center of NGC 2681. 

The geometric properties of the model are fixed by imposing a match with the observed 
isophotes as obtained with the ellipse IRAF task. The photometric profile is derived as a func- 
tion of the free parameters given below. This model is then convolved with the proper PSF — a 
normalized mean of both UV and red-leaked PSFs constructed by means of Tiny Tim (Krist 1992, 
version 4.4) and weighted in accordance to a properly chosen energy distribution — taking correctly 
into account of nonlinearity effects (for FOC images alone; see Appendix of Cappellari et al. 1999). 
The resulting model is then simultaneously fit by means of an iterative procedure to all the original 
images as given by the STScI pipeline (with the exclusion of the highly nonlinear pixels with R < 1" 
in the F342W 1993 image). 

A very good fit is achieved also for the far-UV passbands, (Figure 6), consistent with the lack 
of color gradients. However, only ~50% of the detected flux in the three red-leak-affected FOC 
filters (F175W, F220W and F275W) is genuinely UV light (i.e., falls within the filter FWHM, as 
one can verify by using their transmission curves and our own FOS spectrum), thus weakening 
the importance of this latter outcome, when compared to the evidence quoted above. The perfect 
match obtained with both the 1993 and 1997 FOC observations also assures that the nucleus of 
this galaxy has not significantly varied in surface brightness between the two epochs. However, it 
is worth mentioning that the FOC/F342W profiles of NGC 2681 and NGC 4552 (Figure 7) show 
that the surface brightness of NGC 2681 close to the center is higher by « 4 mag arcsec -2 than 
that of NGC 4552 for R « (X'03, before the variable spike starts dominating its luminosity. As 
a consequence the central UV surface brightness of NGC 2681 is so high that a UV flare of the 
luminosity seen in NGC 4552 would have not been detected with observations of the kind discussed 
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here, leaving the NGC 2681 nucleus apparently quiescent for the observer. 
We find that a unique double power-law (Lauer et al. 1995) 

with parameters 2.0 < a, ft = 1.56 ± 0.02, 7 = 1.19 ± 0.02, R b = l'/l ± Gf.'l can reproduce all of 
the eight profiles (see Figures 5 and 6) and also satisfy the constraint that the FOC dark countrate 
(which is also a free parameter for the FOC frames) remains within plausible ranges (Nota et 
al. 1996). Here 7 measures the steepness of the inner profile (I{R) oc R~ J for R <C Rb), (3 is 
the steepness of the outer profile (I(R) oc R~^ for R » i?t>), a is related to the sharpness of 
the transition and 1^ is the scale factor for luminosity. In the profiles a radius of 0.38 pixel — 
corresponding to the average radius within the central pixel — has been assigned to the central 
value of the surface brightness. Due to the small FWHM and the good sampling of the PSF in the 
UV our FOC/96 non-aberrated UV images provide information down to the central ~ 0.3 pc of 
NGC 2681. 

The profiles in Figures 5 and 6 represent the surface brightness in mag arcsec -2 in the STMAG 
magnitude system. The magnitude zeropoints — computed using the SYNPHOT package within 
IRAF/STSDAS — take into account the 1.27 x higher sensitivity of the FOC zoomed mode, relative 
to the non-zoomed observations (Cappellari et al. 1999; Nota et al. 1996). It is evident from the 
profiles from the 1997 FOC images (Figure 6, right panels) that there is a change of slope of the 
observed profile around R ~ Qf.'l which is not reproduced by the model. This radius corresponds to 
the distance where the nucleus becomes more asymmetric (see Figure 8). This, in turn, suggests that 
some substructure could be present. However, no difference in color or spectral energy distribution 
is observed for this central feature, thus justifying our choice of resorting to a single component 
density profile when modeling the dynamics of the central regions. 
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5. DYNAMICAL CONSTRAINTS ON THE NUCLEAR MASS DISTRIBUTION 

In order to construct a dynamical model of the innermost regions we need to extract a reliable 
intrinsic density profile, starting from the observed photometry. To this aim the analytic double 
power-law photometric profile derived in Section 4 is deprojected by means of an Abel inversion 
(Binney & Tremaine 1987), under the assumption of spherical symmetry. The profile is converted 
to the V-band using the color ^342 — [iy ~ 1.0 mag measured from the nuclear FOS spectrum of 
Figure 1 (the STMAG photometric system has the same zero point as that of the V magnitude band) 
and applying a foreground extinction Ay = 0.07 (Burstein & Heiles 1984; Schlegel, Finkbeiner, & 
Davis 1998). 

The deprojected luminosity density profile in units of L Q y pc -3 is presented in Figure 9. The 
derived density in the central pixel (« O'.'Ol) corresponds to po > 2.4 x 10 6 L Q y pc -3 . This is quite 



-8- 



a high value for the central luminosity density. For comparison, for a sample of 41 spirals analyzed 
by Carollo & Stiavelli (1998), the highest value of the density at 10 pc from the nucleus is a factor 
w 2 lower than the corresponding value we derive for NGC 2681 (pio pc ~ 6 x 10 3 L y pc~ 3 ). 

Velocity dispersion and rotation curve profiles are not available for NGC 2681. Moreover, 
as this galaxy contains three concentric bars (Erwin & Sparke 1999), developing a fully three- 
dimensional model of its mass distribution is problematic, at best. However, experience has shown 
that, while detailed observations and general dynamical models are needed to establish the presence 
of possible BHs in galaxies, a simple isotropic dynamical model gives reasonable order of magnitude 
estimates of the BH parameters. A classic example is the «5x 10 9 M BH in M87 inferred by 
Sargent et al. (1978) using spherical isotropic models, an estimate reasonably close to the most 
recent values (« 3 x 10 9 M ) determined by Harms et al. (1994) and Macchetto et al. (1997) using 
HST spectra of the gaseous disk surrounding the nucleus of M87. 

A detailed analysis of spherical models with cuspy power-law inner surface density profiles 
I(R) oc R-i is available from Dehnen (1993), Tremaine et al. (1994) and Zhao (1996, 1997). The 
dynamics of isotropic power-law cuspy nuclei (without BH) are dramatically different, depending 
on whether 0<7<lorl<7<2(7 being the inner slope of the power-law). For the lower 7 case, 
termed "weak cusp" by Tremaine et al. (1994), the mean-squared line-of-sight velocity dispersion 
approaches zero near the center. For the other case, termed the "strong cusp" regime, the relation 
<Tp oc i? 1 " 7 diverges near the center (finite resolution will actually transform these mathematical 
singularities to finite observable values if 7 < 1.5). 

In either case, if a BH is present oc R v near the center. Thus it is clear that, while with 
< 7 < 1 an increase of in the very inner regions (R <C Rb) is already an indication of a possible 
BH, this is not the case when 1 < 7 < 2, as is such for NGC 2681 (7 ~ 1.19). A BH becomes 
finally undetectable for 7 — > 2, when its presence induces the same trend as a simple increase of 
the M/L stellar ratio. 

Figure 10 shows the kinematical data available for this galaxy: ground based measurements 
of the stellar central velocity dispersions of Dalle Ore et al. (1991, a = 111 ± 16 km s _1 within a 
1'.'5 x 4" aperture) and Schechter (1983, a = 105 ± 9 km s _1 within a 4'.'2 diameter aperture), and 
the gas velocity dispersion of Ho et al. (1997, a = 113 ± 5 km s^ 1 within a 4" x 2" aperture). In 
Figure 10 we also plot our own gas dispersion FOS measurement of a = 204 ± 5 km s" 1 within the 
0'/21 x 0'/21 aperture. The central a is found to increase by a factor of ~ 2 when narrowing the 
aperture from R ss l"5 (ground) to R » O'/l (HST/FOS). 

To investigate this further, we have computed a simple spherical isotropic model for the stellar 
kinematics, starting from the double power-law photometric profile I(R) determined above, and 
assuming a constant stellar mass-to- light ratio T = M Q /L Qv . With these assumptions the pro- 
jected stellar velocity dispersion at a distance R from the center can be predicted from the relation 
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(e. g. Tremaine et al. 1994): 

The stellar density p(r) and the mass within a given radius M(r) have been obtained from the 
deprojected double power-law profile. R is the radial distance on the plane of the sky, while r is 
the radial coordinate measured from the center of the galaxy. In the case a BH lies at the center 
of the galaxy we have to substitute the enclosed mass M(r) in Equation (3) with the expression 
[M(r) + M.], where M. is the mass of the BH. 

In order to compare the predictions with the observations we calculate the velocity dispersion 
<7 a (-R) as measured through a circular aperture of radius R centered on the galaxy nucleus (seeing 
and PSF effects are not modeled here, since they are smaller than the apertures used): 



J R al(R')I(R')R' d Rf 
I(R')R' dR' 



<(R) = u £ ' ' ■ ( 4 ) 



We evaluate this equation for the two cases with and without a central BH. In the former case 
we impose that below R = 0'.'005 (the central resolution of our photometric model, corresponding to 

0. 38 FOC/96 pixels) the stellar density profile flattens, becoming constant (an extreme assumption, 
which maximizes the possible BH mass). We do this in order to avoid the integral at the numerator 
in Equation (4) diverging. To compare the models with the data we further assume that the gas 
behaves as a "hot" system, being in distinct clumps in isotropic motion, with a density proportional 
to the stellar density. This is appropriate if the gas is produced by the stars themselves, a not 
unreasonable assumption. 

Figure 10 shows the results of our modelling with the inclusion of a central BH of mass 
M. = 6 x 10 7 M Q (solid line) and without it (dotted line). It is apparent that a model without 
a BH cannot explain the observations. However, we stress that, due to the above assumption on 
the flattening of the central density profile, the value M. < 6 x 10 7 M Q represents a strict upper 
limit of the mass of a possible central BH in this galaxy. In this respect, one should be aware that 
only measurement of stellar kinematics using long slit STIS spectra will permit us to remove the 
present degeneracy in the modelling of NGC 2681. 

We also note that the ground-based data constrain T to be 0.7 and 1.0 for the two above cases, 

1. e. with and without BH, respectively. A value T « 1 is consistent with simple stellar population 
models of solar metallicity and ages 1-2 Gyr, for both Salpeter and Scalo IMF (Maraston 1998). 

The upper limit for the BH mass determined in this paper (M. < 6 x 10 7 M ) is consistent 
with the recent very tight correlation between the central velocity dispersion a c and the BH mass 
as given by Merritt & Ferrarese (2000). The central dispersion a c « 111 km s" 1 (corrected to r e /8) 
we observe for this galaxy, corresponds in fact to an expected BH mass of M, «8x 10 6 M . 

We remark that in the above BH mass estimate we have modeled the gas as a "hot" system of 
point masses. This appears justified from the fact that within an R l'.'5 aperture the dispersion 
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observed from the gas and the stars is the same (see Figure 10). The gas motion however may 
become more ordered closer to the nucleus. As an extreme case, one could even assume that, at 
HST resolution, the gas lies in a thin disk in circular motion in the combined potential of the BH 
and the stellar potential. 

To test this possibility we have tried to explain the observed emission line profile as produced by 
a thin gaseous disk using the techniques we have developed in Bertola et al. (1998), and subsequently 
integrating the model spectrum within the FOS aperture. Using different BH masses and various 
disk inclinations, we have verified that, under the thin disk hypothesis, it becomes impossible to 
produce an emission line having the width a ~ 204 km s _1 observed with FOS without producing 
a double peaked line profile. Intrinsic dispersion in the gas is thus important and the disk cannot 
be thin. This assures that the adopted spherical assumption for the geometry gives the correct 
order of magnitude estimate for the BH mass. 

6. CONCLUSIONS 

Our analysis of the light distribution in the inner 0.6 kpc of NGC 2681 shows that a double 
power law light profile fits its light distribution both in the NICMOS images and in our set of 
FOC/96 UV images. The good fit of the multiband data obtained with the same light distribution, 
coupled to the close similarity of our FOS spectrum to the IUE spectrum taken with a much larger 
aperture, argues against the presence of sizable color gradients in the center regions of NGC 2681. 
A remarkable consequence of the steep, innermost profile of NGC 2681 is that it translates — when 
observed in the near-UV — into a central surface brightness some 4 mag arcsec~ 2 higher than that 
of the recently identified UV- variable spike of NGC 4552 (Cappellari et al. 1999). This implies 
that a UV-bright flare of comparable luminosity as that observed at the center of NGC 4552 would 
escape detection in the nucleus of NGC 2681. 

The FOS spectrum of the central region reveals the presence of broad emission lines. Compar- 
ison to ground based spectra shows that the gaseous emission is centrally concentrated, most of it 
coming from a region inside the central ~ 6 pc of NGC 2681. On the basis of the observed ratio of 
the emission lines, the nucleus of NGC 2681 falls in the region of LINERs. This suggests a slow, 
steady feeding of a central BH. 

The presence of the BH is more strongly indicated by the increase of a factor ~ 2 of the gas 
velocity dispersion, when narrowing the aperture from R ~ 1'.'5 to R ~ Qf.'l. Such an increase cannot 
be accounted for by the stellar potential alone, as derived by deprojecting the observed surface 
brightness profile. We also find that the gas in the inner regions of NGC 2681 is not supported by 
rotation, as in the case of a Keplerian disk, but mainly by random motions. Assuming that the gas 
can be modeled as a "hot" system, we derive a firm upper limit to the BH mass of M. < 6 x 10 7 
Mq from an isotropic dynamical model. Such an estimate is consistent with the M, — a relation of 
Ferrarese & Merrit (2000) and Gebhardt et al. (2000). 
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In this respect, NGC 2681, with a central velocity dispersion of 111 km s _1 (corrected to r e /8) 
and a predicted BH mass M, ~ 8 x 10 6 M , must be regarded as an object of particular interest 
as it falls in the still poorly populated low-u region of the tight M, — a relation. One should be 
aware, however, that only the knowledge of the radial dependence of the velocity dispersion (and 
not its central value alone, as is our current knowledge for NGC 2681) assures the possibility of 
pinpointing the BH mass in NGC 2681 instead of establishing its upper limit, as we have done here. 
To acquire such data requires HST/STIS slit observations. 

Our dynamical models also indicate a stellar mass to light ratio for the central population 
of 0.7-1, which is appropriate for a 1-2 Gyr old, solar metallicity stellar population. This is fully 
consistent with previous IUE and optical spectroscopy (namely, the measured Mg2 line strength 
index) indicating that the population of this galaxy is dominated by coeval stars formed in a fairly- 
well defined starburst ~ 1 Gyr ago (Burstein et al. 1988). 

Finally we remark that a HI survey of NGC 2681 and its associated galaxies would be highly 
desirable. Such a survey might reveal evidence that NGC 2681 underwent a tidal interaction with 
another galaxy in the recent past, lending support to the hypothesis of a major merging event 
affecting what we see in the center of this galaxy today. On the one hand, NGC 2681 belongs to 
a galaxy group, an environment in which mergers are particularly favored to occur. On the other 
hand, NGC 2681 does not show clear signatures of the recent accretion of a gas-rich dwarf, such 
as morphological distortions or duplicity of the nucleus. While the presence of a uniform, dustless, 
young stellar population in its central 0.6 kpc is puzzling, perhaps we can view NGC 2681 as the 
older brother of M82, seen face-on, after the star burst has run its course. 

MC is grateful to the European Southern Observatory for its kind hospitality in the period 
during which much of the HST data analysis has been completed. DB acknowledges support by 
NASA/STscI through grants GO-03728.01-91A and GO-06309.01-94A. 
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Fig. 1.— Combined G270H, G650L and G780H FOS nuclear spectrum of NGC 2681 within the 
(X'21 x (/.'21 aperture. All spectra have been smoothed with a Savitzky-Golay algorithm (Press et 
al. 1992). The absorption and emission features discussed in the text are marked. 
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Fig. 2. — Detailed comparison between the UV portion of the FOS spectrum within the Of! 21 x Of! 21 
aperture (Savitzky-Golay smoothed) and IUE spectrum of NGC 2681 (Burstein et al. 1988) within 
the 10" x 20" aperture (thick line). The IUE spectrum has been rescaled by multiplying its flux 
by the ratio between the FOS spectrum average in the 5400-5600 A range and the V-flux within 
the IUE aperture quoted by Burstein et al. (1988) for NGC 2681. An additional rebinned FOS 
spectrum (thin line) is superposed to the original spectrum to permit an easier comparison with 
the IUE energy distribution at the same resolution. Note the relative absorption line strengths of 
Mg II (2800 A) and Mg I (2852 A), which indicate that the main sequence stars in this part of the 
galaxy are late A. 
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Fig. 3. — Upper panel: Template-subtracted nuclear FOS G650L spectrum of NGC 2681. The up- 
per, middle and lower plots show the original spectrum, the vertically shifted template constructed 
as discussed in the text, and the difference between the above two spectra, respectively. Only emis- 
sion at H/? is apparent in this relatively low signal-to-noise spectrum. Lower panel: Same as in the 
upper panel for the grating G780H. Superposed to the lower, template-subtracted spectrum is our 
resulting model including the [O I], Ha, [N II] and [S II] emission lines. The individual Gaussian 
emission components of the model are shown at the bottom of the panel. 
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Fig. 4. — Location of the NGC 2681 nucleus (large gray square) on the diagnostic diagrams of Ho et 
al. (1997), as derived from emission line ratios measured from the FOS spectra. The observational 
errors are comparable with the size of the smaller symbols. The large open square marks the 
position of the NGC 4552 nucleus, as determined by Cappellari et al. (1999). The other symbols 
represent the nuclei included in the Ho et al. sample (crosses = H II nuclei, filled squares = Seyfert 
nuclei, filled circles = LINERs, open circles = transition objects). The vertical and horizontal 
lines delineate the boundary adopted by Ho et al. between (a) H II nuclei, (b) Seyfert galaxies, (c) 
LINERs and (d) Transition objects. 
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Fig. 5. — The surface brightness profile of NGC 2681 on the STMAG photometric system, mea- 
sured from the F160W NICMOS image. Diamonds represent the observed values, while the solid 
line represents the double power-law profile of the best fitting photometric model, determined as 
described in the text. Note that exactly the same double power-law profile has been used to model 
the FOC photometry in Figure 6. 
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Fig. 6. — STMAG photometric profiles in the FOC/96 passbands. Dashed lines always represent 
(apart from a scale factor) the unique double power-law model, fitting simultaneously all the HST 
profiles. Such a unique model translates into different expected profiles (solid lines) once the instru- 
mental effects affecting each different FOC image are applied. Diamonds represent the measured 
profiles after subtracting the most appropriate dark level. Error bars represent the allowed range 
of surface brightness at each measured radius, given the possible range of background count rate 
and measured random errors for each image. The fact that in most cases diamonds lie close to the 
upper bar limit comes from the fact that the dark-rate was generally close to the nominal (lowest) 
value. The same double power-law profile also reproduces the NICMOS data in Figure 5. The 
highly nonlinear pixels within m 1" in the F342W 1993 image cannot be reproduced by the model. 
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Fig. 7. — Comparison between the 1996 FOC/96 COSTAR-corrected photometric profiles of 
NGC 4552 (triangles) in the F342W band and NGC 2681 (diamonds) in the F342W+F2ND band 
from our 1997 FOC observation. For each profile (shown in the STMAG photometric system) open 
symbols represent the observed points, the dashed line is the double power-law profile of the best 
fitting photometric model of the galaxy and, finally, the solid line represents the model after the 
inclusion of FOC instrumental effects. A point-like source has been added to the NGC 4552 model 
to reproduce the observations (see Cappellari et al. (1999) for details). 
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Fig. 8. — Contour image of the central arcsec in the 1997 FOC/96 non-aberrated image of 
NGC 2681. This image has been obtained by accurately aligning the three F175W, F275W+F1ND 
and F342W+F2ND FOC/96 images by cross-correlating the frames and then coadding them to 
improve the signal-to-noise. The resulting image has been finally boxcar-smoothed over 3x3 pixels 
((y/04 x Of! 04). The contours are logarithmically spaced. Note the clear asymmetry of the nucleus 
and the irregular appearance of the surrounding regions. 
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Fig. 9. — Deprojected density profile of NGC 2681. The profile has been obtained by Abel inversion 
of the best-fitting analytic double power-law profile. The F342W+F2ND band calibrated profile 
has been first converted to the V band by subtracting the value ^342 — /xy ~ 1.0 mag, measured 
from the FOS spectrum. 
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Fig. 10. — Aperture velocity dispersion (cr a ) profile (dotted line) computed by deprojecting the 
observed double power-law photometric profile and by solving the hydrostatic equilibrium equation 
under the assumption of spherical symmetry and isotropy of the velocity-dispersion tensor. With 
T = 1.0 M /L y this profile agrees with the stellar dispersion measured by Dalle Ore et al. (1991, 
filled square) and by Schechter (1983, filled triangle). The open circle and open square indicate 
the gaseous dispersion measurement by Ho et al. (1997) and our own FOS value, respectively. The 
solid line represents the same isotropic model with the addition of a 6 x 10 7 M central BH and 
T = 0.7. Since a double power-law brightness profile with 7 > 1 and a central BH would generate 
an infinite observed dispersion, in this computation we have assumed that the stellar density profile 
becomes constant inside the FOC/96 resolution limit (see text for details). 
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Table 1. FOS Observation Log. 



Grating 


A range 


Spectral Resolution 


Date 


Exp. 




(A) 


(FWHM A) 




(s) 


G270H 


2222-3277 


1.89 


1997 Feb 2 


2600 


G650L 


3540-7075 


23.4 


1997 Feb 2 


600 


G780H 


6270-8500 


5.26 


1997 Feb 2 


1720 



Table 2. FOC f/96 Observation Log. 



COSTAR 


Observing Mode 


Filter 


FoV 


Scale 
(" pixel- 1 ) 


Date 


Exp. 

(s) 


No 


512zxl024 


F175W 


22"x22" 


0.0225 


1993 Nov 5 


2793 a 


No 


512zxl024 


F220W 


22"x22" 


0.0225 


1993 Nov 5 


1077 


No 


512zxl024 


F275W 


22"x22" 


0.0225 


1993 Nov 5 


597 


No 


512zxl024 


F342W 


22"x22" 


0.0225 


1993 Nov 5 


417 


Yes 


512x512 


F175W 


7"x7" 


0.01435 


1997 Feb 1 


1134 


Yes 


512x512 


F275W+F1ND 


7"x7" 


0.01435 


1997 Feb 1 


546 


Yes 


512x512 


F342W+F2ND 


7"x7" 


0.01435 


1997 Feb 1 


396 



a Sum of two equal exposures 
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Table 3. Line Emission Fluxes and Modeling Parameters 



Spectrum 


Line 


Line Flux 


FWHM 


LOS v 


(1) 


(2) 


(3) 


(4) 


(5) 


G270H 


C II] A2326 


21.9±2.4 a 


770±70 


390±30 


G780H 


[0 I] A6300 


43.9±6.3 


480±ll b 


575.6±5.0 b 


J) 


[N II] A6548 


144.7±2.4 








[N II] A6584 


434.0±7.2 C 






« 


Ha 


118.0±6.4 






J) 


[S II] A6716 


56.1±6.5 








[S II] A6731 


62.2±6.5 







a This line is actually a multiplet. Since the components are unresolved we have modeled it with 
a single Gaussian component. 

b The FWHM and redshift have been constrained to be the same for all lines in the G780H 
grating. 

c The ratio of the two [N II] lines has been forced to be exactly equal to the theoretical value 3.0. 



Note. — Col. (1): FOS grating used in the observations; Col. (2): Identification of the emission 
line; Col. (3): Integrated flux in the emission line in units of 10~ 16 erg s _1 cm ~ 2 ; Col. (4): Measured 
Gaussian FWHM (cry/ 8 In 2) in units of km s -1 ; Col. (5): Observed redshift in km s -1 . 



